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We investigate some fundamental features of a class of non-linear relativistic lagrangian field 
theories with kinetic self-coupling. We focus our attention upon theories admitting static, spherically 
symmetric solutions in three space dimensions which are finite-energy and stable. We determine 
general conditions for the existence and stability of these non-topological soliton solutions. In 
particular, we perform a linear stability analysis that goes beyond the usual Derrick-like criteria. 
On the bcisis of these considerations we obtain a complete characterization of the soliton-supporting 
members of the aforementioned class of non-linear field theories. We then classify the family of 
soliton-supporting theories according to the central and asymptotic behaviors of the soliton field, 
and provide illustrative explicit examples of models belonging to each of the corresponding sub- 
families. In the present work we restrict most of our considerations to one and many-components 
scalar models. We show that in these cases the finite-energy static spherically symmetric solutions 
are stable against charge-preserving perturbations, provided that the vacuum energy of the model 
vanishes and the energy density is positive definite. We also discuss briefly the extension of the 
present approach to models involving other types of fields, but a detailed study of this more general 
scenario will be addressed in a separate publication. 
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The non-linear Born-Infeld (BI) field theory was orig- 
inally introduced to remove the divergence of the elec- 
tron's self-energy in classical electrodynamics Q. The 
BI procedure defines a new lagrangian for the electro- 
magnetic field, as a given function of the two field in- 
variants, which reduces to the Maxwell lagrangian in the 
low-energy limit and exhibits central-field static solutions 
which are finite-energy and stable. This procedure has 
been extended and widely used in many different contexts 
in Theoretical Physics. As examples of these extensions 
we mention: a) The generalization of the BI procedure to 
non-abelian gauge fields and to higher space dimensions, 
suggested by some results in the context of M theory five- 
brane [5| . b) The search for self-gravitating soliton solu- 
tions of Einstein equations 3] . c) The description of dark 
energy in the context of Cosmology, as a gauge field gov- 
erned by a Bl-like action, or as scalar fields with several 
kinds of derivative self-couplings Q . d) The elaboration 
of chiral models with soliton solutions for the phenomeno- 
logical description of the nucleon structure, as an alterna- 
tive to the usual Skyrme model with a stabilizing quartic 
term in the lagrangian . Almost all of these generaliza- 
tions introduce modified lagrangians which are obtained 
from the initial ones using the same BI prescription. Be- 
sides historical reasons, this choice is supported by some 
nice properties of the original and extended BI actions, as 
the existence of finite-energy (electrostatic and dyon-like) 
solutions [6] , the exceptional properties of wave propaga- 
tion (in the electromagnetic original version [7]) or the 
improvement of the stability behavior. It is clear that the 
physical motivations for studying field theories inspired 
on the BI approach are many-fold. Consequently, it is of 
considerable interest to establish under what conditions 



these field theories admit stable particle-like solutions. 
The aim of the present work is to determine these condi- 
tions and to use them to characterize the set of models 
supporting the aforementioned type of solutions. As we 
shall see, the alluded conditions define exhaustively this 
set of models and allow for their explicit determination. 
This provides a wide panoply of supplementary tools in 
order to address the aforementioned problems and others. 

We have first approached these questions in the sim- 
pler case of scalar fields {(j){x^)), with lagrangian den- 
sities defined as arbitrary functions of the kinetic term 
{X = d^(f).d^4)). This is the natural restriction to the 
scalar case of the more general problem outlined for elec- 
tromagnetic fields, where the general lagrangians can 
be defined as arbitrary functions of the field invariants. 
Moreover, the analysis of this problem and its extension 
to the case of many-components scalar fields, are basic 
steps for the later generalization to the case of abelian 
and non-abelian gauge fields. However the results for 
gau ge models will be presented in separate publications 



The search for this kind of scalar models exhibiting 
static soliton solutions in three space dimensions, by cir- 
cumventing the hypothesis of the Derrick theorem 8|, 
has been already partially discussed in an old paper 9|]. 
Here we go beyond these results by establishing the gen- 
eral conditions for the existence of such solutions and per- 
forming a general analysis of their linear stability, beyond 
the Derrick {necessary) criterion. Most of this analysis 
can be easily generalized to other dimensions. 

Let us start with the lagrangian density 
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where the function /(X) is assumed to be continuous 
and derivable in the domain of definition (f2). For the 
purposes of the present analysis we shall call "class-1" 
field theories the models H]) for which f{X) is defined and 
regular everywhere (il = 5ft) and "class-2" field theories 
those with C 5ft, G 17 and open and connected. For 
obvious physical reasons other models are excluded. In 
all cases the associated field equations take the form of 
a local conservation law (in what follows we denote as 
f{X) and f{X) the first and second derivatives of f{X), 
respectively) 



f [X] 



(2) 



The energy density obtained from the canonical energy- 
momentum tensor is 



p^2f(X)[^] -f{X). 



(3) 



We require this energy density to vanish in vacuum and to 
be positive definite everywhere. This imposes the follow- 
ing supplementary restrictions on the Lagrangian den- 
sity: 

/(O) = ; f{X) > (VX) ; f{X) < (VX < 0) 

^(^) > (VX > 0). (4) 

We shall call "admissible" the field models satisfying 
these requirements. For static spherically symmetric 
(SSS) solutions (pir), Eq.^ can be integrated once, lead- 
ing to 



r'cl> = A, 



(5) 



where (f) = d(j)/dr, and A is the integration constant. 
From this equation and the conditions Q we see that 
4> (r) (if unique) must be a monotonic function. Strictly 
speaking Eq.(l5|) determines the field strength for r > 
only. If we substitute the solutions of this equation in 
([2]) we obtain a Dirac 6 distribution of weight 4ttA. We 
can then identify this parameter as the "source point 
charge" associated with the scalar SSS solution (alter- 
natively, if (j) (t) is asymptotically coulombian, A can be 
interpreted as the total scalar charge, continuously dis- 
tributed in space with a density V^0, in analogy with the 
definitions in the non- linear BI models [l[). 

The total energy of the SSS solutions obtained by in- 
tegration of dH]) reads 

poo 

£(A) = -47r / r^f(~(l>'^{r, A))dr = K^'^e{A = 1), 

"'0 

(6) 

where the last equality is a consequence of the invari- 
ance of the solutions of Eq.Q under the scale trans- 
formations 4>{r,t) — > \^^4>{\r^Xt). The convergence of 
this integral depends on the behavior of the SSS field 




FIG. 1: Different possible behaviors of the admissible la- 
grangians with finite-energy static central field solutions. 
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strength at the origin and at r 
power law expression (<f> (r) r'') in both cases 14| . 



it follows from Eq.([5]) that the integrand in ^ behaves 
as r^/(— ^(r)) ~ (/) (r) ~ r'^ and thus the convergence 
of the energy integral requires q > —1 when r ~ and 
(7<— lasr— s-cx). Moreover, the parameter q determines 
the behavior of the lagrangian density f{X) around the 
values of AT = — ^(r) in the limits of the integral. This 
allows to obtain supplementary conditions to be imposed 
on this function in order to have finite-energy SSS solu- 
tions. In this way the discussion of the different possi- 
ble behaviors of the solutions at r ^ and as r ^ oo 
leads to a classification of all admissible models support- 
ing this kind of solutions (see reference ll| for more de- 
tails). When r '--^ we can distinguish three subcases (see 
Fig.l). If > g > —1 (caseA-1) the lagrangian density 
shows a vertical parabolic branch as X — > — oo and the 
field strength of the SSS solution diverges at the origin, 
but the integral of energy converges there and the field 
potential (/)(r) remains finite, li q ~ ( case A-2 ) the field 
strength of the SSS solution remains finite at the origin 
((/) (r ~ 0) ~ C — /3r"'), and the slope of the lagrangian 
density diverges at X = —<j> ^(0) — —C^. The lagrangian 
/(X) itself can take a negative finite value at this point 
(for a > 2; the scalar version of the BI model corresponds 
to cr = 4) or show a vertical asymptote there (for a <2). 
The case q > ( case A-3 ) must be discarded. Indeed, 
if 2 > q > the lagrangian density diverges in vacuum 
and if g > 2 the energy density becomes negative around 
X = 0. 

Concerning the asymptotic behavior of the solutions 
we can also distinguish three subcases (see Fig.l). When 
—2<q<—l ( case B-1 ) the slope of the lagrangian den- 
sity vanishes around X ~ 0. The definiteness of the la- 
grangian around the origin restricts the admissible values 
of q to the rational numbers, by means of the relation q = 
where S ~ N1/N2 is the irreducible ratio of two 
odd naturals such that Ni < N2. For q = —2 (case_^::2) 
the lagrangian density behaves as the D'Alambert la- 
grangian around X '--^ (/(X — > 0^) ~ X) and the field 
strength is asymptotically coulombian. When q < —2 
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( case B-3 ) the slope of the lagrangian diverges at X 0, 
but the lagrangian itself remains well defined there if the 
admissible values of q are restricted to the set of ratio- 
nal numbers through the formula q = SJ^jj where now 
S — iVi/A^2 is the irreducible ratio between an even nat- 
ural A^i and an odd natural N2 > Ni. 

Let us give three illustrative examples of admissible 
models belonging to the different cases and showing soli- 
ton solutions (see Ref. Uj for a more detailed analysis of 
these models): 



Li 
L2 
L3 



odd 3 



X 
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1 



1/2 < a < 1 



(7) 



odd odd 1 

o ',0<a— <l;a — — 7-, > a + - 

2(1 + ^2^)° even odd 2 



The first lagrangian (Li) defines a two-parameter family 
of class- 1 field theories which fall into the cases A-1 and 
B-2 (field strength with an integrable singularity at the 
center and coulombian asymptotic behavior). The con- 
stant A gives the intensity of the self-coupling and the 
family reduces to the D'Alambert model as A ^ 0. It can 
be generalized to larger families of models with soliton so- 
lutions, including lagrangians which arc odd-powers se- 
ries expansions in X or odd functions of X. The second 
lagrangian {L2) defines a two-parameter family of class- 
2 field theories, belonging to the cases A-2 (finite limit 
of f{X) a.t X — ~(p 2(0) and field strength finite at the 
center) and B-2 (coulombian asymptotic behavior). The 
value a = 1/2 in this family corresponds to the scalar 
version of the BI model. As the parameters ^ or 
a ^ 1 all the models converge to the D'Alambert one, 
and the SSS solutions to the Coulomb field. The third la- 
grangian (L3) defines a three-parameter family of class-2 
field theories belonging to the cases A-2 (showing a verti- 
cal asymptote a.t X — — </> ^(0)) and to the cases B-1, B-2 
or B-3, depending on a = 1. In all cases the strength of 

the field is finite at the origin (0 (0) = l//i) and behaves 
like 7--2/(2q!-i) asymptotically. 

Let us now address the question of the stability of the 
SSS finite-energy solutions. The linear stability of these 
solutions requires their energy to be a minimum against 
small perturbations preserving the scalar charge of the 
soliton. At the first order the (vanishing) modification of 
this charge by a small perturbation S(j)(f), assumed to be 
finite and regular everywhere, takes the form 

AA = y d3rV.[/(Xo)(V<5^) - 

-2/(Xo)(W-V<50)V0] =0, (8) 

where now Xq — — (/)2(r). The condition AA = im- 
poses restrictions on the behavior of the admissible per- 
turbations at r = and as r — > 00. In particular, S4> 



must vanish asymptotically faster than the solution (j) 
approaches its asymptotic value. In this way the per- 
turbed fields remain inside the space of functions defined 
by the prescribed boundary conditions which determine 
uniquely the solutions. Introducing the perturbed func- 
tion in the integral of Eq.Q and expanding up to the 
second order, it is easily seen that the first variation of 
the energy vanishes. This is the necessary condition for 
the soliton energy to be an extremum. The second vari- 
ation reads 



dsr /(Ao)(VJ0)2-2/(Ao)(V0.V(50)2 , (9) 



which, owing to Eq.®, converges for any charge- 
preserving perturbation. The positivity of A2S, which 
is the sufficient condition for linear stability, requires 



f{Xo) + 2Ao/(Xo) = 



-2A 
30" (r) 



> 0, 



(10) 



to be satisfied in all the range of values oi X = Xq cov- 
ered by the solution (the equality is obtained by deriv- 
ing Eq.lO with respect to r). This condition is always 
fulfilled by admissible models with finite-energy SSS so- 
lutions, owing to the monotonic character of (j) (f)- We 
conclude that all finite-energy SSS solutions of admissi- 
ble models are linearly stable against charge-preserving 
perturbations. We can perform a more detailed analy- 
sis of the dynamics of the small perturbations {S4>{r,t)) 
governed by the linear equation 



fiXo)V{6(j)) - 2/(Xo)(V0.V(50)V0 



d_ 

di 



fiXo) 



0, 



(11) 



obtained from the linearization of Eq.([2]) around the SSS 
solutions. Note that (fTTj) also takes the form of a con- 
servation law. The general conditions imposed on the 
admissible models, besides the knowledge of the behav- 
ior of the SSS solutions around r = and as r — > cxj, 
allow to perform the standard spectral analysis of these 
linear problems, without explicit specification of the par- 
ticular form of the lagrangians [11| . We summarize here 
the main conclusions of this analysis, which can be de- 
duced from the separation of space and time variables in 
Eg. pT)) together with the boundary condition ([8]). There 
is, in all discrete spectrum of eigenvalues, whose 

associated eigenfunctions are mutually orthogonal and 
finite-norm, with respect to the scalar product defined as 
the spatial integral of the products of the functions with 
f{Xo{r)) as kernel. The corresponding time dependence 
is oscillatory. Any initial perturbation, bounded in the 
norm, remains bounded as time evolves, confirming the 
linear stability. 

All these results can be generalized to the case of iV- 
components scalar fields with a dynamics governed by 
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lagrangian density functions of the form 



N 



(12) 



where, as in the scalar case, f{X) is a given continuous, 
derivable and monotonically increasing function. As al- 
ready mentioned, the analysis of this problem is also a 
necessary step in the generalization of these methods to 
non-abelian gauge field theories. The field equations as- 
sociated with the lagrangians (fT^ take now the form of 
iV local conservation laws: 



(13) 



where X = X^iLi da(l)i-d"(t)i. For the SSS solutions (j)i{r), 



these equations have iV first integrals of the form 



(14) 



where (p^ = d(j)i / dr and A.; are integration constants. The 
canonical energy density reads 



N 



p(x) = 2/(A)^(^) -/(A), 



dt 



(15) 



and is positive definite under the same conditions con- 
straining the function /(A) in the one-component case 
(in fact the set of conditions (H]) defining "admissibil- 
ity" are assumed to hold also in this case). In or- 
der to solve the system (|T4)) we introduce the functions 
Xi{r) = -4>-^{r) and A(r) = % squaring 

and adding Egs. fM)) we are lead to 



-A/(A) = A, 



(16) 



where A = y X]t=i ^1- We see that this equation is for- 
mally identical to Eq.([5|) and, if the function /(A) is the 
same in both cases, there is a one to one correspondence 
between the solutions of the scalar case (j) (r. A) and the 
spheres S\ of radius A in the A^-dimensional A^-space, 
associated with sequences of solutions of the multiscalar 
case of the form 



:(r,A,) = ^,^'(r,A). 



(17) 



The constants A^ can now be identified as "charges" asso- 
ciated with the different components of the SSS field, in 
analogy with the BI and the one-component scalar cases. 
The energy of these solutions, obtained from the inte- 
gration of (fTS]) . is the same as the energy of their scalar 
counterparts, obtained from There is a degeneracy 
in S'a, obviously related to the rotational symmetry in 
the internal space of the lagrangian (fT^ . Moreover, the 



conditions determining the multiscalar field models with 
finite-energy SSS solutions coincide with those already 
discussed in the scalar case. Concerning the conditions 
for stability of the solutions P?|) , the analysis of the one- 
component case can be straightforwardly generalized to 
the present situation ll| . The final conclusion is that the 
multicomponent SSS solitons are linearly stable against 
any perturbation preserving the scalar charges A^, if 
the associated one component solitons are linearly sta- 
ble (note also that transitions in the degeneration sphere 
S'a are blocked by the charge conservation conditions). 

Some additional comments on the general stability of 
these non-topological solitons are in order. Obviously, 
their linear stability does not guarantee the conservation 
(or even a proper definition) of the "soliton identity" in 
presence of strong "external" fields. As is well known, in 
many examples of topological solitons their presence in 
any field configuration can be detected through the ex- 
istence of associated discrete topological charges, which 
are conserved no matter the intensity of external interac- 
tions. In some few cases (always in one-space dimension) 
explicit exact many-soliton solutions have been found, 
allowing the direct analysis of the dynamics of the sys- 
tem in terms of interacting solitons [I4I . But a satisfac- 
tory general analysis of the interaction of solitons with 
strong external fields in three space dimensions is still 
lacking. However, there are some tentative approaches 
to this question which have been developed in the frame- 
work of the Born-Infeld model. We mention the method 
advanced in Ref . [3] , based on the use of the discontinu- 
ity of the field strength at the center of static BI solitons, 
as a marker of the presence and location of the dynamic 
soliton in strong external fields. Since all our soliton solu- 
tions exhibit similar central singularities, this procedure 
could be generalized to the models considered here. But, 
in any case, the permanence of these solitons in strong 
interactions remains a hypothesis which, at best, is com- 
patible with this method. On the other hand, in the case 
of weak external fields linear stability does imply iden- 
tity preservation. Consequently, linear stability is a basic 
condition for the consistency of low- energy calculations 
of the interaction between solitons and weak fields (or be- 
tween distant solitons). The results of such calculations 
may be interpreted in terms of particle-field (or particle- 
particle) force laws and give also a first approach to the 
radiative behavior in these processes 

As already mentioned, the present analysis has been 
generalized to electromagnetic and non-abelian gauge 
fields of compact semisimple Lie groups [lo| . The main 
result amounts to establish a correspondence between 
any scalar model with stable SSS solutions of finite en- 
ergy and families of gauge field models exhibiting similar 
solutions (which can be explicitly written in terms of 
the scalar ones). Conversely, the families so defined 
exhaust the class of gauge field models supporting the 
aforementioned kind of solitons. 
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